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Massive Neutrinos in Cosmology 
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The roles of massive neutrinos in cosmology — in leptogenesis and in the evolution of mass density fluctuations 
— are reviewed. Emphasis is given to the limit on neutrino mass from these considerations. 
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1. INTRODUCTION 

There are two places in cosmology where the 
mass of neutrinos would play a significant role. 
One is leptogenesis in the early universe, and the 
other is the evolution of mass density fluctuations 
that are explored with cosmic microwave back- 
ground (CMB) fluctuations and the formation of 
cosmic structure. There is one more place where 
the presence of neutrinos is generally important 
- primordial nucleosynthesis, but the effect of 
neutrino mass is negligible unless it is unrealisti- 
cally heavy. 

In this talk I discuss at some length the effect of 
the neutrino mass on the CMB fluctuations and 
cosmic structure formation, and the limit on the 
neutrino mass therefrom: we saw a substantial 
progress in our understanding after the reports of 
Wilkinson Microwave Anisotropy Probe (WMAP) 
PP and Sloan Digital Sky Survey (SDSS) [2] but 
also some confusions still remain. We start with 
a brief mention about the status of leptogenesis. 



2. LEPTOGENESIS 

One of the promising ideas for baryogenesis is 
the generation of baryon asymmetry via leptogen- 
esis from the Majorana mass term in the presence 
of the action of sphalerons of electroweak interac- 
tions The simplest scenario assumes delayed 
decay of a thermally produced heavy Majorana 
neutrino to a scalar particle and a light neutrino, 
and its conjugate. For the delayed decay scenario 
to work there is a lower limit on the smallest 
heavy Majorana neutrino mass, as known from 



the argument of GUT baryogenesis j lloKij . which 
in turn leads to an upper limit on the mass of 
light left-handed neutrinos via the sea-saw mech- 
anism. The latest analyses with the Boltzmann 
equation yield |7I8I9| 

m Vi < O.feV, (f) 

for all species of neutrinos. 

It is also desirable to impose the condition that 
pre-existing baryon asymmetry that may have 
arisen from some baryon or lepton number vio- 
lating processes at, e.g., the Planck mass scale, if 
any, is erased so that the prediction of leptogene- 
sis does not depend on the initial condition. This 
demands that the Majorana neutrino be smaller 
than a certain value, so that lepton number vi- 
olation takes place fast enough compared with 
the expansion of the Universe. This gives a lower 
limit on some effective light neutrino mass 



O.OOIeV < m Vm 



(2) 



The combination of the Yukawa terms that ap- 
pear in lepton number violation does not give 
the mass term that is written in terms of light 
neutrino masses; hence the interpretation of this 
lower limit needs some care. This limit does not 
mean that the lightest neutrino must satisfy this, 
but it suggests that the lepton number violation 
process is unlikely to erase pre-existing lepton 
asymmetry if all neutrino masses are as light as 
those violating this limit. 

There is a recent focus of interest if this thermal 
leptogenesis scenario is viable in the world with 
supersymmetry, with which the reheating tem- 
perature cannot be sufficiently high to produce 
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needed heavy Majorana neutrinos without over- 
production of gravitinos. For a recent review, see 

ma. ' 

Experimental verifications of the leptogenesis 
scenario need: 

(1) neutrinos are of the Majorana type; 

(2) find some evidence for the presence of the uni- 
fication scale that is relevant to massive Majorana 
neutrinos (M > 10 9 GeV). More specifically, this 
is the scale where rank of the unifying group steps 
down to 4 at which extra U(l) gauge group is 
broken (see §9.3.3 and §9.6.1 for a detailed 
argument). 

(3) CP violation. However, it would be even more 
a surprise if CP is conserved (i.e., mass matrix is 
real) for some reasons. 

3. COSMIC STRUCTURE FORMATION 
AND MASSIVE NEUTRINOS 

By now we believe we understand the evolution 
of the universe as a whole and of cosmic structure 
at large scales. The latest most important step to 
modern cosmology was the discovery of fluctua- 
tions in the cosmic microwave background by the 
COBE satellite in 1991. This indicated that we 
are on the right track to understand the cosmic 
structure formation. At the same time it gave 
compelling evidence for cold dark matter (i.e., 
dark matter that was non-relativistic when it was 
decoupled from the thermal bath) that dominates 
the matter component of the Universe. In the 
last ten years the cosmological paradigm was also 
shifted. The cosmological constant was an anath- 
ema in 1990, but now all observations converged 
to pointing towards the existence of the cosmo- 
logical constant that dominates the energy of the 
Universe. It is surprising that no observations 
have yielded evidence against the A dominated 
cold dark matter (ACDM) universe. The re- 
cent results from WMAP and SDSS strongly sup- 
ported this standard view based on the ACDM 
Universe, as demonstrated by the convergence of 
the cosmological parameters {12113114] to 

if =71±5km s^Mpc" 1 , 

O m = 0.28 ± 0.03, n A = 0.72 ± 0.03, 

Q m + A = 1.01 ±0.01 . (3) 
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Figure 1. Observational power spectrum from 
CMB fluctuations, galaxy clustering, gravita- 
tional lensing shear, cluster abundances and Ly- 
man a clouds, compared with the prediction of 
ACDM model. The figure is taken from the SDSS 
paper, Tegmark et al. ^5] 



We also know that baryons amount to only 1/6 of 
O m , and only 6% of them are comprised in stars 
and stellar remnants. 

Another important measure is the power spec- 
trum that characterises matter fluctuations, 

P(k) = J d 3 xe tk - x (6{x)5(Q)), (4) 

where S(x) is the density contrast at the position 
x. A variety of observations at vastly different 
cosmological epochs (CMB at z = 1090, galaxy 
clustering at z ~ 0.1 — 0.4, gravitational lensing at 
z ~ 0.5 — 1, cluster abundances at z ~ 0.1 — 0.2), 
when scaled to z — 0, yield P(k) that is described 
by |fc|"T(fc) with n = 1 ± 0.02 and the trans- 
fer function T(k) predicted by the ACDM model 
with the cosmological parameters specified above 
|15I16| (see Figure 1). This consistent description 
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of P(k) constitutes additional evidence that sup- 
ports the standard model, nearly scale-invariant 
initial adiabatic perturbations (S = 3ST/T) grow- 
ing by gravitational instability in the ACDM uni- 
verse. We must add a remark that no models are 
known, other than inflation, that generate these 
fluctuations consistent with observations, while 
a successful model of inflation from the particle 
physics point of view is yet to be found. 

The massive neutrinos contribute to the mass 
density of the Universe by the amount 

,_ 2 



n„ = 



(5) 



94.1 eV 

where h — 0.71 is the conventional notation for 
the Hubble constant. In the following discussion 
we assume three species of light neutrinos and do 
not consider exotic possibilities which are occa- 
sionally discussed by particle physicists. The neu- 
trino mass and its limits discussed in what follows 
are of the order of 1 eV, which is much larger than 
than 0.05 eV derived for the mass difference for 
the limiting case of the hierarchical-mass neutri- 
nos. It is therefore appropriate to consider three 
degenerate neutrinos, and we assume this in our 
considerations. 

3.1. Effects of massive neutrinos on the 
evolution of cosmic fluctuations 

A successful model is obtained for cosmic struc- 
ture formation without having massive neutrinos. 
This means that massive neutrinos only disturb 
the agreement between theory and observations, 
hence leading to a limit on the neutrino mass. 

The well-known effect of massive neutrinos is 
relativistic free streaming that damps fluctua- 
tions within the horizon scale. One electron volt 
neutrinos are still relativistic at matter-radiation 
equality, which takes place atTwl eV, and then 
tend to smear fluctuations up to ^100 Mpc co- 
moving scale, thus diminishing the power of P(k) 
for these scales. This effect becomes stronger as 
the cosmological mass density of neutrinos, hence 
the neutrino mass, increases (see Figure 2) (see 
|17I18| ). Therefore, the empirical knowledge of 
P(k) across large to small scales gives a constraint 
on the summed mass of neutrinos. 

The best determination of P(k) for large scales 
is given by CMB multipoles observed by WMAP 
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Figure 2. Power spectrum predicted in the 
ACDM model including massive neutrinos with a 
degenerate mass specified in the legend. The top 
curve represents the case of massless neutrinos 
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PP. There are a number of ways to obtain P(k) 
at small scales: the use of galaxy clustering, the 
cluster abundance, the Lyman a cloud absorption 
optical depth, and the cosmic shear due to gravi- 
tational lensing. Limits on the neutrino mass are 
derived from a single or combined use of these 
observations. One may also input other observa- 
tions, which serve to restrict various cosmologi- 
cal parameters, which may indirectly improve the 
limit on the neutrino mass in a multidimensional 
parameter search. 

(1) Limits derived from galaxy clustering varies 
from £nv < 0.6 eV [E] to 2.1 eV [E] at a 
95% confidence, both authors using 2dFGRS data 
(see also [201 ) 1 - The limit using the SDSS data 
is X) m f < I-? e V EH- (See also |23 , where 
J2 fn v < 0.75 eV is concluded using SDSS and 
2dFGRS.) While I quoted the limits all obtained 
by using only CMB and galaxy clustering here, 
the origin of different upper limits is not clear. It 
is possible that the use of SDSS or 2dFGRS would 
leads to different limits, as the fall of P(k) to- 
wards smaller scales is more gentle with the 2dF- 
GRS, thus giving a stronger constraint on m„. 
A suspect is the convergence of practical applica- 
tions of the Markov chain Monte Carlo method to 
calculate the likelihood, which we discuss briefly 
in section EP1 

(2) The use of the cluster abundance enables us to 
estimate the mass fluctuations that are the quan- 
tity relevant to us. This was made in the early 
work |22j . but no updates were attempted after 
WMAP. 

(3) With the Lyman a cloud absorption power 
derived from fluctuating optical depths, 2 one can 
explore P(k) at the smallest scale (23124) . giv- 
ing the strongest limit m v < 0.42 eV (95%) 
|25| (see [Zjo\ for the earlier work). This result, 
however, is more model dependent in the sense 



1 Allen et al. I2UI concluded a finite value for the neutrino 
mass at 2<r. Here we take their upper 2a value as a limit, 
which is 1 eV. 

2 Lyman a clouds are the objects that cause Lyman a ab- 
sorption lines with the neutral hydrogen column density 
10 13 — 10 17 cm -2 in the quasar spectrum when they in- 
tervene the line of sight to a quasar. The 'clouds' are 
identified with moderate overdensities of the hydrogen gas 
of temperature ~ 10 4 K, governed by photoionisation and 
adiabatic cooling. 



that one must invoke simulations to extract P(k) 
from observed flux power spectrum, which suf- 
fer from significant uncertainties associated with 
modelling and simulations. 

(4) With gravitational lensing one can directly 
explore mass fluctuations. For the moment the 
statistical error is not sufficiently small, but this 
provides us with a promising method. This may 
also be used to set the normalisation of the power 
spectrum derived from galaxy clustering. 

Elgar0y and Lahav [2T give a summary of lim- 
its obtained in the literature. 

We note that a decrease of spectral index n 
(red tilting) has an effect similar to that of mas- 
sive neutrinos for small scales. Hence, unless the 
length scale explored spans a wide range or n is 
constrained, m„ and n are positively correlated, 
i.e., an 'artificial' blue tilting would be cancelled 
by an appropriate value of m v ; this weakens the 
constraint on m v . Elgar0yetal. |28| derived from 
2dFGRS power spectrum 'alone' m v < 2 eV, 
but this is under the assumption that the spec- 
trum is exactly flat, n = 1. Allowing for n > 1, 
a sensible limit is lost. We definitely need CMB 
data to derive a limit on the neutrino mass. 

3.2. Cosmological limits on the neutrino 
mass: caveats 

The methods listed above all suffer from dif- 
ferent systematic errors. What we really need 
is the mass fluctuation £ m = {5 m (x)6 m (0)) in 
the linear regime. We expect that nonlinearity 
for the length scale relevant to galaxy clustering 
> 10 Mpc is only modest (e.g. [2H]), but with the 
use of galaxy clustering we must assume that the 
galaxy distribution traces mass allowing for some 
constant factor, i.e., £ s (x) = b 2 t; m (x), where b is 
left usually as a free parameter, representing "bi- 
asing" of galaxies relative to mass. This is proba- 
bly not a bad approximation at the accuracy that 
concerns us at the present, but for higher accu- 
racy we do not know how good is this approxi- 
mation: at least we know that biasing depends 
on morphology and luminosity of galaxies (e.g., 
POj). Hence the sample used must be homoge- 
neous in this regard. The systematic difference 
in P(k) between the two groups must also be re- 
solved. It should also be noted that errors for 
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the galaxy clustering data given in the literature 
may not represent properly the systematics that 
are associated with the evaluation of the window 
function of the survey data analysis. 

The use of the cluster abundance is proposed 
as a way to estimate the mass fluctuation tak- 
ing advantage that the cluster mass can be esti- 
mated with various means, from velocity disper- 
sion measurements, X ray data, and gravitational 
lensing. For the moment, however, there are non- 
negligible uncertainties in the cluster mass esti- 
mate. The disadvantage of this method is that 
one needs a very large cluster sample to estimate 
P(k) for varying scales. 

One can explore mass fluctuations in the small- 
est scale with the use of Lyman a cloud absorp- 
tion power spectrum; thus one can get potentially 
the information most sensitive to m„. The prob- 
lem is that one needs substantial corrections to 
unfold the matter power spectrum P(k), which 
can be done only with simulations. While a great 
success of numerical study was to enable us un- 
derstand the nature of Lyman a clouds |31| (see 
|32| for a review), it is not easy at a quantitative 
level to document the systematic errors. The sim- 
ulation always suffers from mesh effects, and this 
is particularly true when baryons are included. 
With the inclusion of baryons one has to deal 
with stars and their feedback, which is a difficult 
astrophysical problem. Specifically for the case 
of Lyman a clouds, we must be concerned with 
uncertainties as to the heating rate from ionis- 
ing radiation and thermal history of clouds. The 
estimate of the mean level of absorption optical 
depths is an additional source of errors. While 
much work is currently being done, the assess- 
ment of errors in P(fc) derived from Lyman a 
clouds is not easy a problem. 

In principle the use of cosmic shear is the best 
for estimating the mass fluctuation. The nonlin- 
earity correction is still needed, but it can be done 
with simulations without baryons, which are rel- 
atively simple. A difficulty is that the signal is 
small compared with noise so that one needs a 
very large galaxy sample to attain good statis- 
tics, and that one needs to know instrumental 
distortions of images accurately. 
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Figure 3. CMB multipoles for the MDM like so- 
lution suggested by Elgar0y & Lahav (J2 m v = 4 
eV, dotted curve) are compared with those for the 
best ACDM solution (with massless neutrinos, 
solid curves). The % 2 is higher by 50 with this 
MDM like solution. The data are from WMAP. 



3.3. Limits on neutrino mass from CMB 
alone? 

There is a controversy as to whether a sensible 
limit is derived on the neutrino mass from CMB 
multipoles alone. This consideration is meaning- 
ful since one might be able to derive the most 
robust limit on the neutrino mass free from sys- 
tematics that are difficult to control with data on 
galaxies. Tegmark et al. |13| . however, showed 
that no limits are obtained from CMB alone: 

m v < 9 eV is allowed at one sigma. Elgar0y & 
Lahav 133] endorsed this result, emphasising the 
need to employ large-scale structure data. They 
also suggest that there is a solution that repro- 
duces CMB multiples only with CDM and mas- 
sive neutrinos that satisfy Q m + 0„ = 1 without 
a cosmological constant, i.e., the 'old' mixed dark 
matter (MDM) scenario is viable (see Figure 3), 
with the only price being a low Hubble constant. 

On the other hand, Ichikawa, Fukugita & 
Kawasaki [3J] claimed that a sensible limit 

m v < 2 eV (at a 95% confidence) can be 
derived from CMB (WMAP data) alone. They 
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argued that the effects of massive neutrinos on 
CMB multipoles cannot be absorbed into the shift 
of cosmological parameters (notably the Hubble 
constant) if neutrinos are already nonrelativistic 
at the recombination epoch, i.e., m v > 3T roc ~ 
0.55 eV. The principle is that the presence of 
non-relativistic neutrinos causes a decay of grav- 
itational potential in pre-recombination epoch, 
which amplifies the acoustic oscillation that ap- 
pears as the second and third peaks (while the 
first peak receives little the effect from the poten- 
tial decay, since the multipoles corresponding to 
free streaming of the m„ > 0.55 eV neutrino are 
I > 250 which are larger than the position of the 
first peak, l x = 220). 

The shifts of the position and the height of the 
first peak, dominantly caused by modifications of 
the angular diameter distance and of the inte- 
grated Sachs- Wolfe effect, respectively, are still 
significant for m„ < 0.55 eV, but the heights of 
the second and third peaks relative to the first 
are hardly modified. These effects can be ab- 
sorbed into the shift of the other cosmological 
parameters. When the second and third peaks 
change in addition for m v > 0.55 eV, however, 
one cannot absorb the effects into the shift of 
cosmological parameters. This is the basic mech- 
anism how the constraint on m v is derived from 
CMB multipoles. The Hubble constant prior as 
claimed in is not essential. The numerical 
analysis of Ichikawa et al. showed that the in- 
crease of x 2 with the neutrino mass is very slow 
till m u ~ 1-5 cV and gave m n < 2.0 eV at 
a 95% C.L. 

Figure 4 compares the CMB multipoles for the 
case of m Vi = 2.0 eV with the best model 
having massless neutrinos that gives the global 
X 2 minimum. This demonstrates the accuracy 
of the WMAP data, which has a power to dis- 
tinguish the two curves. It is argued that 
one cannot obtain a limit significantly better than 
0.55 x 3 ~ 1.6 eV even if the precision of CMB 
multipoles is improved. 

Note that an increase of the peaks at higher I 
modes with massive neutrinos mimics blue tilt- 
ing, and hence there emerges a negative correla- 
tion between n and m v , which is opposite to that 
expected from the effect on the power spectrum, 
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Figure 4. CMB multipoles for the ACDM solu- 
tion with massive neutrinos that is allowed just at 
a 95% confidence (J2 m v = 2.0 eV, dashed curve) 
are compared with those for the best ACDM so- 
lution with massless neutrinos (solid curve). The 
data are from WMAP. 



as explained in the previous subsection. 

Elgar0y & Lahav |23 ascribed the discrepancy 
between Tegmark et al./ Elgar0y & Lahav and 
Ichikawa et al. to the difference of priors, but 
they did not give an analysis that supports this 
statement. I would ascribe the origin to insuffi- 
cient samplings of the Markov-chain Monte Carlo 
methods the former authors adopted to estimate 
the likelihood functions, especially away from the 
point at m„ = 0. When a chain is not sufficiently 
long the Markov-chain Monte Carlo methods of- 
ten result in likelihood functions that differ in 
shapes depending on the initial condition. This 
was the reason why Ichikawa et al. preferred to 
use a deterministic algorithm to search for the 
minima at given m„ with successive parabolic ap- 
proximations, with the further use of the Vegas 
integration code to check for the likelihood func- 
tion. In this approach, priors do not play an im- 
portant role. They also remarked that they could 
not reproduce the likelihood function of Tegmark 
et al. as a function of n away from the global min- 
imum, indicating that their likelihood function 
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may probably be incorrectly evaluated, whereas 
the behaviour around the global minimum was 
reproduced very well. 

Ichikawa et al. further searched for MDM like 
solutions, but failed to find any with A% 2 smaller 
than 16. The curves shown in Figure 3 is the 
solution suggested by Elgar0y & Lahav which 
gives x 2 higher than the best solution by A% 2 = 
50. Although the shape of the multipoles appears 
similar to that of the best ACDM solution to eyes, 
the accuracy of the WMAP data rules out the 
MDM model at a high confidence. 

3.4. Conclusions and prospects 

I discussed that there are mainly two effects 
that lead to constraints on the neutrino mass, 
damping of the small scale power due to neutrino 
free streaming, and decay of the gravitational po- 
tential in the pre-recombination epoch that am- 
plifies the acoustic oscillation in high multipole 
modes. 

The limit on the neutrino mass from cosmic 
fluctuations is £ m Vi < 2 eV (95%) from CMB 
alone, which is subject to systematic errors to 
the least extent. The assumption is that ACDM 
model is correct and the power spectrum obeys 
a power law, allowing for a small departure from 
the power law to the extent as expected in slow- 
roll inflation. The limit obtained by combin- 
ing the galaxy clustering with the CMB data is 
J2 m vi ^5 1 — 1.7 eV (95%). Unfortunately, the er- 
ror arising from the Monte Carlo sampling is not 
well documented, and an accurate limit is yet to 
be found for a given set of the clustering data. A 
systematic error in the power spectrum between 
two major groups is also yet to be resolved. The 
current most stringent limit is < 0.42 eV 

(95%) derived with the use of Lyman a cloud ab- 
sorption power spectrum by unfolding the matter 
power spectrum with the aid of simulations of the 
Lyman a cloud. It is not easy to assess its relia- 
bility at the moment. 

The accuracy of the CMB multipole data will 
be greatly increased, especially at the higher mul- 
tipoles with the launch of Planck [HSj- However, 
a straightforward analysis of its data is unlikely 
to bring a drastic improvement in the constraint 
on the neutrino mass. Ichikawa et al. considered 



that the limit that can be achieved from CMB 
alone will be J2 m ^i < 1-2 eV. Efforts to remove 
systematic errors are needed for small scale clus- 
tering data. Particularly important is to enhance 
the reliability of simulations in a way convincing 
to everybody and to document errors involved in 
the results of simulations. 

Clearly, the goal is to give a constraint as 
small as X) TO ^ ~ 0.05 eV, which is the mini- 
mum mass indicated by neutrino oscillation ex- 
periments. Kaplinghat et al. -37 proposed to 
use gravitational lensing signals in CMB polari- 
sation. They forecast that one would get a limit 
^2 Tn-ut < 0.15 eV with Planck, but need a new 
project particularly sensitive to CMB polarisa- 
tion to get to < 0.04 eV. An accurate simulation 
is also needed to unfold non-linear effects of mass 
clustering in small scales. Wang et al. [HHj pro- 
posed to use large cluster surveys. Their forecast 
is m vi < 0.03 eV by combining cluster sur- 
vey using LSST with multipoles from Planck. 
Here the estimate of cluster mass is crucial, which 
they hope to perform using simulations of clus- 
ters. 
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